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Abstract. The Riesz s-energy of an iV-point configuration in the Euclidean space M. p is defined 
as the sum of reciprocal s-powers of all mutual distances in this system. In the limit s — > 
C ~ ~) , the Riesz s-potcntial l/r s (r the Euclidean distance) governing the point interaction is replaced 

f\J ■ with the logarithmic potential log(l/r). In particular, we present a conjecture for the leading 

term of the asymptotic expansion of the optimal L2-discrepancy with respect to spherical caps 
' on the unit sphere in R ti + 1 which follows from Stolarsky's invariance principle [Proc. Amer. 

Math. Soc. 41 (1973)] and the fundamental conjecture for the first two terms of the asymptotic 
expansion of the optimal Riesz s-energy of N points as N — > oo. 
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1. The discrete Riesz s-energy problem 



43 

The Riesz s-energy (real s ^ 0) and the logarithmic energy (s = 0, by convention) of an iV-point 
^ r~| configuration Xjy with points Xi, . . . , xjy in the Euclidean space W (p> 1) are defined as 

C3 1 N N N N 

^-ggj^, Eo( *„) = gg> g _. 

j^k j^k 

The logarithmic energy can be understood as the limiting case s — > 0. For large s the nearest 
neighbor's interaction dominates and the other limiting case s — > oo gives the best-packing problem 
(or Tammes problem [43]). Most deeply studied is the classical (Newtonian) case s = p — 2 for the 
' harmonic potential l/r p ~ 2 , naturally with an abundance of literature in physics in the Coulomb 

CO . case (p = 3), where the Riesz s-energy essentially is the potential energy of an ensemble of 

' (positive) unit point charges placed at the points of the configuration X^. 

For this note of interest are sums of distances (that is s = —1) for points on the sphere § d 
because of the close connection to the spherical cap L 2 -discrepancy defined in Definition [3] by 
means of Stolarsky's invariance principle (Proposition 0]) and to the worst-case error for equal 
weight quadrature formula for functions in the unit ball in a certain reproducing kernel Hilbert 
space over S d (cf. [10] and [H]). 

The N -point s-energy of an infinite compact set A C R p is defined as follows: 

£ s (A;N):=sup{E s (x u ...,x N ) : xi,...,xjv <£ A} if s < 0, 

£ s (A;N):=M{E s (xi,...,x N ) : xi,...,xjv £ A} if s > 0. 

It gives the optimal s-energy an iV-point configuration Xn in A can assume; that is 

£ S (A; N) = E S (X^) for an optimal s-energy system xff. 

Let d denote the Hausdorff dimension of A. When normalized appropriately (such that the total 
'charge' of the iV-point configuration is one) the quantity S S (A\N)/N 2 remains finite (and has a 
limit as N — > oojj in the potential-theoretical regime s < d but grows beyond any bound in the 
hypersingular case s > d. In the regime s < d tools and result from classical potential theory 
(see, for example, Bjorck [5] (s < 0), Saff and Totik [37] (logarithmic case) and Landkof [5S]) are 
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'Indeed, for example, for s > the sequence N(N — 1)/£ S (A; N) is monotone and bounded and converges to 
the 'generalized transfinitc diameter of A of order s' defined by Polya and Szego |35j . 

1 



2 



J. S. BRAUCHART 



used to treat the limit N — > oo. The limit distribution of optimal s-energy configurations is given 
by the (unique if — 2 < s < d) equilibrium (or extremal) measure ha maximizing (if s < 0) or 
minimizing (if s > 0) the energy integral 



^a\tA-= J J Mx,y)dMx)dMy), 



where fc s (x,y):=|x — y\~ s if s ^ and fco(x,y):= — log |x — y|, among all Borel probability 
measures fi supported on A. The extremal quantity V S (A) = Xj/i^] is called the s-energy of the 
set A. If s > d, then T s \p] = +oo for any measure [i supported on A, that is V S (A) = +oo. In 
this case methods and results from geometrical measure theory are successfully used to gain more 
insight. One outcome is that other normalizations of the optimal s-cncrgy would give a converging 
sequence for certain classes of compact sets (see, for example, Kuijlaars and Saff [28] and Hardin 
and Saff [25]). These results are the building blocks for the fundamental conjecture below. 

The discrete Riesz energy problem is concerned with investigating properties of sequences of 
optimal (and nearly optimal) Riesz s-energy configurations. Questions concern (i) explicit compu- 
tations of optimal configurations, (ii) limit distribution, (Hi) asymptotic expansion of the optimal 
energy (cf. Smale's Problem #7 [40]), (iv) geometric properties ('well-separation', 'mesh'-norm). 
We refer the interested reader to the survey articles Saff and Kuijlaars [36] and Hardin and Saff [24] . 

Let S d denote the unit sphere in The following long-standing open fundamental conjecture 

for the asymptotic expansion of the optimal Riesz s-energy is known. (In the upcoming paper [T^j 
this conjecture and its history will be discussed in more detail.) Let Hd denote the d-dimensional 
Hausdorff measure (normalized so that the d-dimensional unit cube in W has measure 1). 

Conjecture 1 (Fundamental Conjecture). Let d>2 and —2 < s < d+ 2 (s ^ 0,^0 Then 

£ s (S d ; N) = V s (5 d ) N 2 + s ' d N 1+s / d + 7^(S d ; N), 

[H d (S d )] s/d 

where lZ s (E> d ; N)/N 1+s / d ~ e — > as N — > oo for some e > possibly depending on d and s. 

This conjecture combines known results on the basis of the principle of analytic continuation: 

• In the potential-theoretical regime — 2 < s < d (s 7^ 0) the dominant term grows like N 2 
and the coefficient is given by the s-energy of the unit sphere 

For other values of s it is understood to be the analytic continuation of the right-hand 
side above to the complex s-plane except at the simple poles at s = d, d + 2, . . . , 2d — 2 
(finitely many) if d is even and at s = d, d + 2, d + 4, . . . (infinitely many) if d is odd. 

• In the hyper-singular case s > d the leading term behaves like N 1+s / d and it is shown 
by Kuijlaars and Saff [35] that the limit £ s (S d ; N)/N 1+S ' d existsQ It is believed that 
the constant C s ^i can be analytically continued to the s-plane, which is supported by the 
conjecture for its value in case of d = 2 also provided in [28] . 

• In the boundary case s = d [28] gives that 

^(S d ;iV)~^SiV 2 logiV asJV^, 
rid\P ) 

which can be also understood by means of a limit process s — > d assuming that the poles 
of V s (E> d ) and C s .d at s = d cancel each other. (Obviously, one can also consider limit 
processes s — > s' , where s' is a pole of V s . In such a case one would gain information about 
the singularity at s = s' of the coefficient of the corresponding term in the asymptotic 
expansion of the optimal s-energy (s near s') provided this lesser-order term exists. The 
limit process s — ¥ connects the asymptotic expansion of the optimal logarithmic energy 
and optimal s-energy for s near 0.) 



^or s < — 2 the extremal distribution is concentrated in no more than d + 1 points on S d ([6]). 
'The existence of this limit is proven for the class of rectifiablc a!-dimcnsional manifolds in |25 | and for weighted 
Riesz s-energy for rectifiable d-dimensional sets in [7]. 
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The constant H d (M d ) /H d (E> d ) is the ratio of the volume of the unit ball M d in M. d and the surface 
area of S d (denoted by u> d ) and can be expressed in terms of a ratio of surface areas or a ratio of 
gamma functions 

(2) n d (S d )-d oj d ~ d ^Y{d/2) (asd^oo). 
It is a long-standing open problem what the precise value of the constant 

(3) C s . d = lim £ a ([0,l]*;N)/N 1+ '' d 

is for d > 2. (In the case d = 1 Martfnez-Finkclshtcin et al. showed that G s \ = 2£(s), where £(s) 
is the Riemann zeta function.) The significance of the constant C s ^ d follows from the fact that 
the limit lim s _ i . 00 [Cs^] 1 / 8 can be expressed in terms of the 'largest sphere packing density' in R d , 
which is only known in three cases: d = 1, 2 and d — 3 (Kepler conjecture, rather recently proved 
by Hales [23] )• For d = 2 Kuijlaars and Saff [28] obtained the estimate 

f, s (S 2 ;A) (V3/2) S/2 C A ( S ) 
limsup — — t— < , s > 2, 

and they conjecture that equality holds above; that is: 

Conjecture 2 (Kuijlaars and Saff [28]). For s>2, C s , 2 = (T^)^ 2 f A (s). 

The function Ca( s ) i s the zeta function of the hexagonal lattice A = {m (1, 0) + n(l/2, \fZ/2) : 
m,n € Z}, which was used in [28] to locally approximate a minimal s-energy configuration on the 
sphere in the hypersingular case s > 2. For Res > 2 it is defined as 

Ca( s ) = |7p = X! ; iW2 

and it is understood that Ca( s ) is the mcromorphic extension to C of the right-hand side above. 
The function Ca( s ) admits a factorization (cf., for example, jTS] Chapter X, Section 7]) 

(4) Ca(s)=6C(s/2)L„ 3 (s/2), Re S >2, 

into a product of the Riemann zeta function £ and the first negative primitive Dirichlet L-Series 

(5) L_ 3 ( S ):=1-1 + 1-1 + 1-..., Rc,s>l. 

Interestingly, it is assumed in physics that a hexagonal configuration has lowest potential energy 
but there seems to be no mathematical proof for this. However, it is known that the hexagonal 
lattice is (even universal) optimal among all lattices in M 2 , see Montgomery [35] and Cohn and 
Kumar [TB] , 

2. Uniform Distribution and Discrepancy 
Definition 1. A sequence {Ajv} is asymptotically uniformly distributed on E> d if 

,. # {k : x fc ,w £ B} 

for every Ud-measurable clopen set B in § d . 

Informally speaking: a reasonable test set gets a fair share of points as N becomes large. 
The quality of a sequence {AT at} of Appoint systems can be quantified using the discrepancy 

#{fc:x fe eB} 



D{JF; xi, . . . , xa?):= sup 
Be/ 



A 



<Td(B) 



measuring the maximum deviation between the uniform measure (limit distribution of optimal 
configurations) and empirical point distribution with respect to a family T of test sets (for example 
spherical caps). The spherical cap discrepancy of an A-point configuration Xn will be denoted by 
Dc(AOv). 
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A different approach to uniform distribution makes use of numerical integration and function 
spaces. Let 

N 



I[/]:= f fda d , Q w [/] : =I^/( Xfe 

Jsd k=l 



Definition 2. {Xn} is equi- distributed with respect to every function in the space of continuous 
functions if[l 

lim Q N [f] = I[f] for every / e C(S d ). 

N—Hxi 

As it is well-known both Definitions Q] and [2] are equivalent. 

On the sphere one has no satisfactory analogue to the celebrated Koksma-Hlawka inequality 
in the unit cube. Cui and Freeden [T7] introduced the concept of generalized discrepancy on § 2 
based on pseudo-differential operators. For a particular choice D they obtained a Koksma-Hlawka 
like inequality 

|Qiv[/]-I[/]| < V6D CF (X N )\\f\\ a3/2(s2) , 

where / is from a certain Sobolev space H 3 / 2 (§ 2 ) whose reproducing kernel is defined using D. In 
this context, a sequence {Xjy} of TV-point systems is called D-equidistributed with respect to all 
functions in H 3 / 2 (S 2 ) if lim.jv-i.oo ^cf(Xn) = 0. Moreover, the generalized discrepancy associated 
with D has a closed form expressible in terms of elementary functions 

f N 

47T[D CF (X N )} 2 = 1-^j £ log(l + | Xi -x ft |/2) 2 . 

j,k=l 

(Sloan and Womersley [39] showed that [Dcf(Xn)] 2 has a natural interpretation as the worst-case 
error for Qn for function from the unit ball in H 3 / 2 (S 2 ) provided with a norm which is equivalent 
to the one used by Cui and Freeden.) This approach is followed further in [14] leading to the 
generalized discrepancy 

N 



(6) [D(A W )] 2 = i - -L £ 



3 N 2^^-^\ 

k,e=i 

associated with the Sobolev space H 3 / 2 (S 2 ) with the reproducing kernel K(x, y) = (8/3) — |x — y|. 

3. The spherical cap L 2 -discrepancy on E d 

A spherical cap on § d centered at x £ S d is the set 

C(x;<):={ye§ d : (x,y) ><}, 

where (•,•) denotes the usual inner product. The family of all spherical caps on E> d forms a 
discrepancy system (see [20] ) giving rise to the notion of asymptotic uniform distribution 

Definition 3. The spherical cap ^-discrepancy of an TV-point configuration X^ on S rf is given 
by 



Dj?(A w ): 



1 Js d 



dcr £i (x) dt 



N 

The optimal ^-discrepancy of iV-point configurations on S d is denoted by 

D^ 2 (S d ; A^):=inf {d^ 2 (Aat) : X N C § d | . 
An obvious upper bound of the L2-discrcpancy is in terms of the spherical cap discrepancy: 

D^(Xjv) < V2 B C (X N ), X N C S d . 
The next result connects the sum of distances, the L2-discrepancy and the distance integral. 



That is, the discrete probability measure associated with tends to as N — > oo (in the weak-* limit). 



OPTIMAL DISCRETE RIESZ ENERGY AND DISCREPANCY 



Proposition 4 (Stolarsky's invariance principle [H]). Let d>2. The 



N 



x fe 



H d {^ d ) 
H d {M d ) 



D^(X 



|x-x'|d<rd(x)d<r d (x'), X N C S d 



Clearly, reduction of the L2-discrcpancy means increase of the sum of distances and vice versa. 
Stolarsky's invariance principle provides a simple way to compute the L2-discrcpancy of a given 
iV-point configuration on S d . In the language of the discrete energy problem we have 

A' 



(7) 



N i 



V-^ d ) 



1 

N~ 2 



x fe 



[^ 1 (X JV )-F_ 1 (S d )iV 2 l/7V 2 



£ to 

j,k=l 

Evidently, the square of the optimal L,2-disrepancy is closely related to the second term of the 
asymptotic expansion of the optimal Riesz (— l)-energy, cf. Conjecture [T] This connection between 
the spherical cap (L2-)discrepancy and the second term in the asymptotic expansion of the s-energy 
is exploited in [SJ, Grabner and Damelin [18]. 

The IL,2-discrepancy of an A-point configuration Xn is minimial if and only if the sum of 
distances of Xn is maximal. Thus, optimal L,2-discrcpancy configurations arc, in fact, maximal 
sum-of-distance configurations or, cquivalently, maximal Riesz (— l)-energy configurations. Such 
configurations can be generated using numerical optimization which is a highly non-linear process. 

Stolarsky used his invariance principle and discrepancy results of Schmidt [38] on the discrep- 
ancy of spherical caps to estimate the difference 

|x-X |dcr d (x)dcr ( i(x ) — 

and obtained the correct order of A for the upper bound. Harman [26] improved the lower bound 
(for the general invariance principle) and, finally, Beck [3] obtained the correct order of A for 
the Euclidean metric. By means of the invariance principle these bound for the energy difference 
translate into bounds for the spherical cap IL,2-discrepancy. 

Proposition 5. Let d > 2. There exist constants c', C > such that for sufficiently large N 
c' AT-V2-i/(2d) < D ^(S d ; N) < C 7V-i/2-i/(2d)_ 

Remark. By Proposition [5] the correct order of the decay of the optimal ]L,2-discrepancy on § d in 
terms of powers of the number of points N is iV -1 / 2-1 ^ 2 ^ which is the same rate (apart from the 
-y/log N term in the upper bound) as for the optimal spherical discrepancy on § d . 

Proposition[5] rises the question if D^, 2 (E> d ; N) N 1 / 2+1 / ( 2d ) has a limit as N — > oo. By employing 
the connection to the discrete Riesz s-energy problem, where one has conjectures regarding the 
asymptotic expansion of the optimal s-energy, we derive the following conjectures for the leading 
term in the asymptotic expansion of the optimal IL,2-discrepancy. 

Conjecture 3. Let d>2. Lf Conjecture^ holds, then 
(8) D^ 2 (S d ;A)~ ^ d Ar-V2-i/(2d) + ... asN->oo, 



-C. 



i. d 



Justification. By Stolarsky's invariance principle (Prop. |4] and Eq. 0) 



n d (s d ) 



D^ 2 (S> d ; N) 



H d (M d ) 

Application of the fundamental Conjecture [T] yields 



n d (s d ) 



n d 

The conjecture follows. 



Dc 2 (§ d ; AO 



[n d mr l/d 



£-i(E> d ; A) — T^_i(§ d ) A 2 ] /N 2 . 

N- 1 - 1 / d -TZ- 1 (E d ;N)/N 2 . 



-C-i.d 



□ 



'The given version is a special case of the result in 1411 which is extended in [8], 
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For the unit sphere § 2 one can make the conjecture more precise. 
Conjecture 4. If Conjectures^ and\^hold, then 

D^ 2 (§ 2 ; N) ~ A 2 ^-V2-i/0M) + . . . as jv ->■ oo, 

where 



A 2 = y - [- C(-l/2)] L_ 3 (-l/2) = 0.44679728350408 . 

Justification. Since Hd(^> d ) = cl>2 = 47r and (by (J2J)) 



H d (S rf ) 



1 r(3/2) _ i 

2 V^r(i) ~ I' 



d=2 

one obtains (using Conjecture [2] and relations (j4j and (|8])) 



A 2 = J- 



1/2 -6C(-l/2)L_ 3 (-l/2) 




(4-)- 1/2 

0.44679728350408 . 



□ 



The case d = 1. Stolarsky' invariance principle also holds for <i = 1 as one can see from the proof 
given in [10]. In this case one has (cf. Eq.s ([lj and ([2])) 

N 



± |x J --x fe | + i[D^) 



2 4 , 
= -, X N C§\ 

7T 



On the other hand (in joint work with Hardin and Saff) |13| we obtained a complete asymptotic 
expansion of the Riesz s-energy C s (N) of the iVth roots of unity which represent optimal iV-point 
configurations for s > —2. In general, for s e C with s ^ 0,1,3,... 

C S (N) = V„(S) N 2 + ^V +s + t^- V a n (s) ((s - 2n) N 1+s - 2n + (N- 1+Be '-^ ). 

C 271 ") ( 27r ) ^ 

The coefficients a n (s), n > 0, satisfy the generating function relation 



smnz 



= J2 a n(s)z 2n , \ Z \ <1, SeC. 



n=0 

It follows that 

p-i 



1 



7T 



Dq 2 2 = 2[ (2 ^ ( ) _ 1 1)] + j^zt E C("l - 2n) iV- 2 - 2 " + G(iV- 2 -^) 



n=l 



where (Bo = 1, -Bi = —1/2, . . . are the so-called Bernoulli numbers and (— l)" +1 i?2n > 0) 
a„(-l) C(-l - 2n) = ( ^ - ^ 2 " < 0, n>l. 

In the last step it was used that the Riemann zeta function at negative integers can be expressed 
in terms of Bernoulli numbers. 

Numeric Results. Recently, in joint work with Josef Dick we investigated the properties of so- 
called digital nets (see Niederreiter [M] and Dick and Pillichshammer [TO]) lifted to the sphere 
by means of an area preserving map. Numerical experiments for Sobol' Sequences lifted to § 2 in 
[TT] suggest that the order of N is optimal but the constant seems to be too large (in the range 
[0.50, 0.59] for N = 2 m , m = 1, . . . , 20). (Optimal (— l)-energy configurations are more difficult to 
come by and the number of points for which numerical optimization is feasible and the result could 
be trusted to be the global maximum is very limited. Cf. 0|4j[5] an d references cited therein.) 
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4. Estimating the spherical cap discrepancy 



The classical Erdds-Turdn type inequality (cf. Grabner |22| . also cf. Li and Vaaler |31| ) estimate 
the spherical cap discrepancy in terms of Weyl sums. Recently, by generalizing LeVeques result 
for the unit circle |30j . Narcowich et al. |33| obtained LeVeque type inequalities on the sphere: 



Cl 



Z{d,l) 

E 

= 1 m=l 



E 



N-l 

-Y 



li, m (Xj) 



1/2 



oo Z{d,i) 

E^E 

f=l m=l 



AT-1 

-E 

AT ^ 



Y iim (xj) 



l/{d+2) 



where af.~T(t~ 1/2) /T{i + d+l/2) x l/£ d+1 =:6^ for some positive constant c\ and C2 and Z(d, n) 
denotes the number of linearly independent real spherical harmonics Yi^ n of degree £. In Sun and 
Chen |42j 'spherical basis functions' (as a counter part to radial basis function on spheres utilized 
in [33]) are used to investigate uniform distribution on spheres. 

G. Wagner obtained ground-breaking results concerning estimates of the Riesz s-energy and 
explored connections between s-energy and discrepancy (see [UJ 25J 25] ) . In fact, the proofs of 
both the Erdos-Turan type inequality and the LeVeque type inequality can be modified yielding 
estimates of the spherical cap discrepancy in terms of s-encrgy. A different approach exploits 
the connection between the error of numerical integration for polynomials (cf. Damclin and 
Grabner [18]) and, by utilizing a result of Andricvskii, Blatt and Gotz [1] (see this author [9]). 

We close this note by recalling that optimal Riesz s-energy configurations may not have optimal 
spherical cap discrepancy (essentially of order 0(N~ 1 / 2 ~ 1 ^ 2d ^), see Beck [3]) which is reflected 
in Korevaar's conjecture |27| claiming for the harmonic case s = d — 1 that the spherical cap 
discrepancy is of order 0(N~ x / d ). This was essentially proved by Gotz [21]. He also showed that 
the bound is sharp for d = 2. The conjecture is still open for other values of s and d. 



Acknowledgement: The author is grateful to the School of Mathematics and Statistics at 
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